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EXPERIMENTAL  CALIBRATION  OF  A  VECTOR  ELECTRIC  FIELD 
METER  MEASUREMENT  SYSTEM  ON  AN  AIRCRAFT 

I .  INTRODUCTION 

The  measurement  of  the  electrostatic  field  in  the  atmosphere  can  be  made 
both  at  ground  level  and  at  altitude.  Ground  based  measurements  can  only 
provide  guesses  as  to  the  field  conditions  occurring  within  a  cloud.  They 
are  not  in  situ  measurements.  Altitude  measurements  can  be  made  with  both 
balloons  and  aircraft.  Aircraft  have  several  advantages  over  balloons. 
Statistical  averages  can  be  performed  over  an  area;  thus  helping  to  filter 
effects  such  as  local  production  of  aerosols  and  space  charges.  Balloon 
measurements,  such  as  performed  by  Hilhleisen  [1971]  are  more  sensitive  to 
these  local  variations.  The  aircraft,  a  more  mobile  and  steerable  platform, 
can  more  easily  investigate  many  geographical  conditions.  The  disadvantages 
of  the  aircraft  are  the  higher  self  charge  and  exhaust  space  charges  caused 
by  the  engines  and  the  difficulty  in  obtaining  the  form  factors  needed  for 
absolute  calibration.  The  costs  also  tend  to  be  higher.  An  aircraft  based 
platform  was  chosen  for  our  measurements. 

Electrostatic  fields  have  been  measured  aboard  aircraft  with  shutter  type 
field  mills  for  several  decades.  Calibration  includes  the  determination  of 
the  enhancement  coefficients  which  relate  the  fields  seen  at  a  particular 
location  on  the  aircraft  skin  to  the  true  external  fields.  Some  of  the 
investigators  are  Gunn  et.  al.  [19A6],  Clark  [1958],  Anderson  [1966],  Imiani- 
tov  and  Chubarlna  [1967],  Fitzgerald  and  Byers  [1958],  Pelton  et.  al.  [1953], 
Bly  and  Nanevicz  [1977],  Christian  et.  al.  [1983],  and  Larouche  et.  al. 
[1985a].  Gunn  et.  al.  [19A6]  and  Imianitov  and  Chubarina  [1967]  did  not 
install  the  required  number  of  A  mills  necessary  in  order  to  resolve  the 
field  into  its  components.  Clark  [1958]  and  Anderson  [1966]  had  four  field 
mill  installations,  but  only  calibrated  two  of  the  three  vector  components. 
Fitzgerald  and  Byers  [1962]  indicate  the  capability  of  calculating  all  three 
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components  of  the  vector  electrostatic  field  and  aircraft  self  charge,  but 
further  scrutiny  indicates  that  several  questionable  symmetry  assumptions 
have  been  made  (e.  g.  -  the  belly  mills  located  near  the  nose  and  the  tail 
have  the  same  sensitivity  to  a  field  directed  along  the  nose-to-tail  axis). 
Also,  one  coefficient  was  not  calculated,  the  authors  stating  that  "further 
flight  measurements  will  be  made."  Furthermore,  no  estimate  of  the  error  of 
the  measured  fields  is  given.  Larouche  et.  al.  [1985b]  and  Bly  and  Nanevicz 
[1977]  also  indicate  resolving  of  the  field  components;  however,  both  cali¬ 
brations  are  based  on  scale  model  measurements  done  in  an  applied  external 
field.  These  measurements  are  suspect  due  to  the  inability  to  accurately 
model  the  exact  field  meter  location  and  exposure.  Pelton  et.  al.  [1953]  use 
sophisticated  equipment  in  order  to  keep  the  self  charge  of  the  aircraft  near 
zero.  They  make  a  reasonable  calculation  of  the  aircraft  self  charge,  but  a 
determination  of  the  enhancement  coefficients  is  apparently  not  performed. 
Jones  et.  al.  [1985]  indicates  that  a  full  vector  calibration  is  near  comple¬ 
tion  for  the  New  Mexico  Tech  aircraft. 

The  statistically  observed  higher  strike  rate  to  aircraft  in  flight  leads 
to  the  hypothesis  that  most  lightning  strikes  to  aircraft  are  triggered  by 
the  intrusion  of  the  aircraft  into  a  region  of  sufficient  electrical  energy 
density.  During  the  summers  of  1984  and  1985,  a  joint  lightning  strike  pro¬ 
ject  [USAF  Flight  Dynamics  Lab  (USAF/FDL) ,  Federal  Aviation  Administration 
Atlantic  City  Technical  Center  (FAA/ACT),  NASA  Kennedy  Space  Flight  Center 
(NASA/KSC),  Office  National  d’Etudes  et  de  Recherches  Aerospatiales  (ONERA) , 
and  Naval  Research  Lab  (NRL)],  was  performed  where  an  instrumented  CV-580 
turboprop  aircraft  was  flown  into  thunderstorms.  NRL’s  installation  con¬ 
sisted  of  four  shutter  field  mills.  This  would  allow  investigation  of  the 
triggering  hypothesis  if  a  full  vector  calibration  is  obtained.  This  report 
focuses  on  the  problems  associated  with  the  calibration  procedure.  A  similar 
less  detailed  report  has  been  given  by  Anderson  and  Bailey  [1986a].  This 
report  supersedes  the  previous  one. 

Because  of  the  considerable  time  and  expense  involved  in  installing  field 
mills  in  an  aircraft,  design  considerations  along  with  the  appropriate  place¬ 
ment  of  the  field  meters  on  the  aircraft  must  have  prime  consideration. 

After  installation,  the  form  factors  must  be  obtained.  The  problems  encount¬ 
ered  during  the  calibration  procedure  will  be  addressed.  A  matrix  of  en¬ 
hancement  coefficients  will  be  defined  in  terms  of  the  field  mill  locations. 


The  determination  of  these  enhancement  coefficients  with  aircraft  maneuvers 
in  a  fair  weather  field  will  be  described.  Error  estimates  are  made  for  all 
the  enhancement  coefficients.  ‘Since  the  ultimate  goal  is  studying  the  fields 
at  the  mill,  the  true  external  fields,  and  the  self  charge  in  relation  to  the 
triggering  hypothesis j  the  enhancement  matrix  must  be  inverted.  A  Monte 
Carlo  technique  will  be  used  to  estimate  the  standard  deviations  associated 
with  the  inverted  matrix,  the  external  fields,  and  the  self  charge.  Finally, 
some  lnq>rovements  that  are  suggested  by  the  analysis  will  be  enumerated. 


II.  FIELD  MILL  DESIGN 

The  calibration  of  field  mills  on  an  instrumented  aircraft  could  be  done 
both  experimentally  and  theoretically.  Estimation  of  errors  are  difficult  in 
theoretical  calculations;  therefore,  an  experimental  calibration  was  per¬ 
formed  for  the  CV-580  aircraft.  Experimental  calibrations  require  a  known, 
constant,  and  homogeneous  field  against  which  to  calibrate.  The  vertical 
fair  weather  field  can  have  these  attributes  along  with  the  added  advantage 
of  only  one  component  to  measure;  therefore,  it  was  chosen  as  our  calibration 
standard. 

The  proper  choice  of  design  characteristics  of  field  mills  can  be  crit¬ 
ical  to  the  calibration  of  an  aircraft  since  fair  weather  fields  can  be  as 
low  as  tens  of  volts  per  meter.  Although  aluminum  is  an  attractive  construc¬ 
tion  material  from  weight  considerations  it  can  have  severe  problems  arising 
from  contact  potentials  if  used  in  the  sensing  areas.  Laboratory  tests  per¬ 
formed  at  NRL  [Willett  and  Bailey,  1983]  indicate  that  non-magnetic  304 
stainless  steel  gives  good  results  with  respect  to  reduced  contact  poten¬ 
tials.  Earlier  experiences  indicate  that  spacings  between  the  rotor  and 
stator  of  several  mm  can  further  assist  in  reducing  the  effects  of  contact 
potentials  without  seriously  reducing  the  sensitivity  of  the  field  mill.  The 
desired  measurement  capabilities  dictate  the  area  of  the  sensor,  the  number 
of  blades,  and  the  rotational  speed  since  these  in  turn  determine  the  field 
meter's  sensitivity  and  frequency  response.  A  wide  dynamic  range  is  required 
if  both  fair  and  foul  weather  fields  are  to  be  measured.  Exposure  of  the 
sensing  area  to  undesired  charges  must  be  eliminated  since  they  will  be 
indistinguishable  from  true  fields;  therefore,  insulators,  which  inevitably 
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build  up  surface  charges,  must  be  guarded  from  the  moving  parts  of  the  sen¬ 
sing  area.  Rotor  grounding  is  also  of  importance  because  synchronous  cur¬ 
rents  flowing  to  the  rotor  through  noisy  or  high  resistance  contact  can  pro¬ 
duce  error  signals.  For  low  level  signals,  a  synchronous  demodulator  is 
necessary  in  order  to  recover  the  signal  from  the  noise.  The  power,  signal 
and  reference  paths  must  be  diligently  shielded  to  minimize  extraneous  sig¬ 
nals.  Finally,  the  power  frequency  should  be  filtered  from  the  signal  and 
reference  paths. 

The  schematic  field  meter  design  and  a  picture  of  the  final  product  are 
shown  in  fig  1.  Dimensions  of  the  CV-580  aircraft  constrained  the  diameter 
to  3.25  inches  (8.26  cm).  A  two  bladed  design  with  a  rotation  speed  of  8000 
rpm  yielded  a  signal  frequency  of  266.7  Hz.  Rotor  grounding  was  provided  by 
silver  graphite  brushes  against  a  coin  silver  slip  ring.  Previous  experience 
indicated  that  brass  and  carbide  brushes  are  inferior  with  respect  to  rota¬ 
tional  noise.  [The  group  at  New  Mexico  Tech  immerses  the  back  of  the  rotor 
shaft  in  a  pool  of  mercury  with  good  success  (J.  J.  Jones,  personal  communic¬ 
ation)].  A  stainless  steel  shell  surrounding  the  sensing  parts  during  opera¬ 
tion  insures  that  contact  potentials  between  the  sensing  parts  and  the  air¬ 
craft  body  are  kept  at  a  minimum.  The  relative  phase  of  the  rotor  with  res¬ 
pect  to  the  stator  is  generated  at  the  back  of  the  motor  unit  by  a  magnetic 
rotor  and  can  be  mechanically  adjusted,  by  the  phase  adjusting  ring,  through 
360  electrical  degrees.  The  overall  design  of  this  field  meter  system  facil¬ 
itates  cleaning  of  the  rotor,  stator,  and  insulators  which  is  essential  to 
operation  in  fair  weather  calibrations.  Time  and  monetary  constraints  prohi¬ 
bited  a  dual-stator  differential  design  that  could  have  doubled  the  signal 
strength. 

The  signals  from  the  field  meter  heads  were  transmitted  along  separate 
anti-microphonic  cables  (RG-149/U  -  essential  in  order  to  eliminate  noise)  to 
amplifiers  inside  the  aircraft.  In  order  to  facilitate  repairs  and  to  con¬ 
serve  space  at  the  field  meter  head,  the  preamplifier  (fig.  2)  was  located  at 
the  amplifier  end  of  the  anti-microphonic  cable.  Since  a  field  mill  is  a 
current  generator,  the  preamplifier  was  configured  as  a  current  sensor.  This 
maintains  the  stator  voltage  near  ground  which  helps  eliminate  leakages  due 
to  voltage  drops  across  the  insulators.  The  preamplifier  is  completely 
shielded  from  the  rest  of  the  amplifier  stages,  is  heavily  decoupled  at  the 
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Fig*  1  -  Detailed  drawing  and  photograph  of  field  meter  head. 
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power  leads,  and  has  back  to  back  diodes  across  Its  Inputs  to  protect  against 
surges.  A  DC  voltage  (hereafter  called  the  bucking  voltage)  can  be  applied 
to  the  stator  through  this  module  to  create  an  artificial  field  at  the  field 
meter  heads  for  sensitivity  calibration. 

The  rest  of  the  amplifier  (fig.  3)  consisted  of  identical  low  pass  fil¬ 
ters  in  both  the  signal  and  reference  paths,  eliminating  both  the  power  freq¬ 
uency  and  signal  harmonic  energy.  The  Frequency  Devices  lowpass  cauer- 
elliptic  filters  had  a  cutoff  of  311  Hz  chosen  such  that  a  null  point 
occurred  at  the  power  fundamental  frequency  of  400  Hz.  The  signal  was  ampli¬ 
fied  and  the  reference  path  was  phase  shifted  by  a  180  degree  variable  half 
lattice  phase  shifter.  This  insured  coherent  phases  of  the  signal  and  refer¬ 
ence  paths  which  were  then  fed  into  an  Analog  Devices  synchronous  demodulator 
(AD630KN).  The  demodulator  has  capabilities  of  narrow  band  transmission  thus 
eliminating  low  frequency  noise  components  and  is  followed  by  single  pole 
smoothing  filters  for  either  thunderstorm  (T  =  0.01  sec)  or  calibration  pur¬ 
poses  (r  =  0.157  sec).  During  thunderstorm  flights  higher  time  resolution 
was  desired  while  during  calibration  the  noise  level  was  required  to  be  mini¬ 
mal.  The  signal  then  passed  through  one  final  stage  of  DC  amplif ication 
before  output  to  the  recording  system. 

III.  RECORDING  SYSTEMS 

The  recording  systems  consisted  of  a  six  channel  Astro-med  (W602XL)  strip 
chart  (frequency  response  =  DC  to  100  Hz)  and  a  twenty  channel  Kyowa  (RTP- 
602A)  FM  analog  recorder  (frequency  response  -  DC  to  313  Hz,  S/N  -  46  dB). 

The  four  field  mill  outputs  were  each  recorded  on  two  ranges  separated  by  a 
factor  of  10.  A  ten  second  slow  time  code  provided  by  USAF/FDL  and  roll  and 
pitch  from  the  aircraft  inertial  navigation  system  were  also  recorded.  The 
measurement  of  roll  and  pitch  are  critical  for  an  accurate  calibration  since, 
especially  during  pitch,  an  accurate  angle  was  difficult  to  hold.  Subse¬ 
quently  these  signals  were  digitized  at  35  and  3.5  samples/sec  for  thunder¬ 
storm  and  calibration  data  respectively. 


IV.  LOCATION  AND  INSTALLATION  OF  FIELD  MILLS 


Since  three  field  components  end  the  aircraft  self  charge  must  be  meas¬ 
ured,  &  minimum  of  four  independent  field  mills  are  required.  Given  the 
coordinate  system  defined  in  fig.  4,  then  a  nodal  curve  can  be  defined  as  a 
curve  on  the  aircraft  skin  where  one  of  the  field  component  contributions  is 
zero  and  a  nodal  point  as  the  intersection  of  two  nodal  curves.  Optimum 
locations  would  then  be  at  nodal  points.  In  principle,  nodal  curves  could  be 
obtained  by  exposing  the  aircraft  to  a  field  in  a  particular  direction  and 
finding  the  locus  of  points  on  the  aircraft  skin  where  there  would  be  a  null 
response.  There  are  three  nodal  curves  and  at  least  six  nodal  points  on  any 
aircraft.  On  the  CV-580  the  nodal  points  occur  near  the  wingtips,  at  the 
nose  and  the  tail,  and  on  the  top  and  bottom  of  the  fuselage  (fig.  5).  The 
simplest  case  is  the  Ey  nodal  curve  (no  response  to  an  Ey  field),  which  is  a 
locus  of  points  along  the  top  and  bottom  of  the  fuselage  created  by  cutting 
it  with  the  one  plane  of  symmetry  perpendicular  to  the  wingspan  and  parallel 
to  the  fuselage.  The  other  nodal  curves  and  points  are  indicated  approxi¬ 
mately  in  the  figure.  No  attempt  has  been  made  to  calculate  their  exact 
locations.  Mazur  et.  al.  [1986]  have  determined,  with  model  calculations, 
approximate  nodal  points  for  the  NASA  F-106  aircraft.  The  curvature  of  the 
aircraft  makes  the  exact  locations  difficult  to  determine. 

A  tail  location  on  the  rudder  is  impractical  and  a  nose  location  may  be 
precluded  because  of  particle  impaction  on  the  exposed  sensor  parts  and  the 
need  for  a  more  powerful  motor  to  overcome  the  increased  airflow  drag.  The 
wingtips  provide  good  measures  of  both  the  wing  to  wing  electric  field  and 
the  aircraft  self  charge.  Wingtip  locations  provide  a  large  field  enhance¬ 
ment  factor  which  makes  measurement  of  the  fair  weather  field  for  calibration 
purposes  much  easier.  Model  calculations  have  larger  errors  associated  with 
modelling  the  wingtip  compared  to  fuselage  locations  since  grid  sizes  are  not 
in  general  fine  enough.  For  model  calculations,  reasonable  locations  might 
be  inboard  on  the  fuselage  if  the  shielding  properties  of  the  wings  are  con¬ 
sidered.  For  those  who  want  to  compare  model  and  experimental  calibration 
techniques,  it  might  be  appropriate  to  have  redundancy  along  the  wingtip  axis 
(two  on  the  wingtips  and  two  inboard).  Finally,  installation  sites  should 
avoid  as  much  as  possible  protrusions  such  as  antennas,  lights,  and  wheels 
which  could  cause  shielding  effects  or  accumulate  surface  charges. 


indicated 


The  field  meter  locations  on  the  CV-580  (fig.  6)  were  two  mills  symmetri¬ 
cally  opposed  on  the  wingtips  (approximately  at  nodal  points)  and  two  mills 
on  the  belly  of  the  fuselage  (along  a  nodal  curve),  as  far  forward  and  aft, 
respectively,  as  possible  in  order  to  maximize  the  separation  for  the  nose- 
to-tall  field.  The  nose  location  was  discarded  both  from  the  necessity  of 
the  radar  Installation  and  from  the  water  Impaction  expected  during  cloud 
penetrations.  Unfortunately,  a  field  meter  at  the  top  of  the  fuselage  was 
not  used  because  of  the  fear  that  it  would  fill  with  water.  Subsequent  ana¬ 
lysis  indicates  that  a  field  meter  at  this  location  would  be  valuable.  In  a 
typical  installation  (fig.  7),  the  field  meters  were  attached  by  four  stain¬ 
less  steel  screws  which  facilitated  their  removal.  In  this  wingtip  configur¬ 
ation  the  glass  light  cover  was  sprayed  with  a  conductive  spray  (Chemtronics 
Static  Free)  during  calibration  flights  to  eliminate  static  charges. 

V.  PROBLEMS  ASSOCIATED  WITH  FAIR  WEATHER  CALIBRATION 

Previous  experiences  at  NRL  on  the  PB4-Y2,  EC 12 IK,  and  the  S2D  aircraft 
indicate  that  with  reciprocating  engines  the  aircraft  self  charge  was  rela¬ 
tively  low.  Boulay  and  Laroche  [1982]  conclude  for  a  Gloster  Meteor  turbojet 
aircraft  little  engine  charging  occurs.  The  CV-580  turboprop,  however,  was 
found  to  exhibit  severe  charging  problems  such  that  the  fair  weather  field 
was  completely  masked.  In  order  to  reduce  the  self  charge  to  a  fairly  stable 
low  level,  .003  inch  (.0762  mm)  stainless  steel  wires  [Anderson  +  Bailey, 
1986b]  were  attached  to  the  trailing  edges  of  the  ailerons  (fig.  8).  This 
reduced  the  aircraft  potential  to  the  corona  potential  of  the  wicks,  low 
enough  such  that  a  sensitive  range  was  reached  by  offsetting  the  remaining 
self  charge  with  a  reasonable  bucking  voltage. 

VI.  CALIBRATION  MATRIX  THEORY 

In  order  to  discuss  the  details  of  the  calibration  procedure,  a  few  defi¬ 
nitions  must  be  made.  Given  the  coordinate  system  (fig.  4),  where  the  air¬ 
craft  coordinates  (dashed)  are  related  to  an  absolute  coordinate  system 
(solid),  a  system  of  equations  describing  the  response  of  the  four  field  mill 
system  can  be  written  as: 


DISCHARGE 

WICKS 


DISC 

WI 


Attachment  of  stainless  steel  wires  to  the  trailing  edges  of  the 
ailerons.  Both  the  standard  discharge  wicks  and  the  attached  wi 
are  shown. 


n  x 


(1) 


Ep  *  Px  Ex  +  Py  Ey  +  Pz  Ez  +  Pq  Q 

Es  *  Sx  Ex  +  Sy  Ey  +  Sz  Ez  +  Sq  Q 

Ef  *  Fx  Ex  +  Fy  Ey  +  Fz  Ez  +  Fq  Q 

Et  =  Tx  Ex  +  Ty  Ey  +  Tz  Ez  +  Tq  Q 

where  Ex,  Ey,  Ez,  and  Q  are  the  field  components  and  the  self  charge  with 
respect  to  the  aircraft  coordinates;  Ep,  Es,  Ef  and  Et  are  the  fields  at  the 
mill  for  the  PORT,  STBD,  FORWARD  (under  the  fuselage),  and  TAIL  (also  under 
the  fuselage)  mills  respectively;  and  the  coefficients  Px  ...  Tq  are  the 
enhancement  coefficients  (note  that  they  are  dimensionless  except  for  the 
self  charge  coefficients).  Symmetry  considerations  suggest  that  the  wingtip 
mills  will  respond  to  an  Ey  field  with  equal  magnitude  and  opposite  sign  and 
to  the  other  field  directions  and  to  the  self  charge  with  equal  magnitude  and 
same  sign.  Responses  to  Ex  and  Ez  fields  are  included  here  because  it  is  not 
yet  known  how  close  we  are  to  wingtip  nodal  points.  Furthermore,  the  belly 
mills,  located  on  the  centerline  of  the  aircraft,  should  be  insensitive  to  an 
Ey  field.  Here  we  are  quite  certain  that  we  are  on  a  nodal  curve  since  this 
is  the  only  true  plane  of  symmetry.  The  matrix  equations  thus  reduce  to: 

Ep  «  Px  Ex  +  Py  Ey  +  Pz  Ez  +  Pq  Q 

Es  =  Px  Ex  -  Py  Ey  +  Pz  Ez  +  Pq  Q 

(2) 

Ef  =  Fx  Ex  +  Fz  Ez  +  Fq  Q 

Et  =  Tx  Ex  +  Tz  Ez  +  Tq  Q 

In  order  to  take  advantage  of  the  wingtip  symmetry,  it  is  convenient  to 
derive  and  record  the  two  quantities: 

Y  =  (Ep  -  Es)/2  =  Py  Ey  (3) 

Q'  =  (Ep  +  Es)/2  =  Px  Ex  +  Pz  Ez  +  Pq  Q  (4) 

which  give  a  direct  measure  of  both  Ey  and  self  charge  (it  turns  out  that  Q' 
is  only  slightly  contaminated  by  the  Ex  and  Ez  fields).  In  the  enhancement 
matrix  above  it  is  assumed  that  the  coefficients  are  statistically  indepen¬ 
dent  and  the  vectors  are  linearly  independent. 
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The  external  fields  and  the  self  charge  are  obtained  by  inverting  the 
matrix.  The  results  are  the  following* 


’FqTz 

TqFz]  fEp  +  Es' 

fTqPz 

FqP  z 

-  - - - -  - 

►  ( -  -  Tz  Ef  + 

•  Fz  -  - 

.  Pq 

Pq  2 

l  Pq 

Pq  . 

COMMON  DENOMINATOR 


(Ep  -  Es) 
2  Py 


'TqFx  FqTx  fEp  +  Es’  TqPx  'FqPx 

.  Pq  Pq  l  2  Pq  J  Pq 


Fx  Et 


COMMON  DENOMINATOR 


(FzTx  -  FxTz) 


’Ep  + 
.  2 


+  (PxTz  -  PzTx)Ef  +  (PzFx  -  PxFz)Et 


Pq  *  COMMON  DENOMINATOR 


where  COMMON  DENOMINATOR  is  given  by 

FqPxTz  TqPxFz  TqPzFx  FqPzTx 


+  FzTx  -  FxTz 


Note  that  the  field  components  depend  only  on  ratios  of  the  self  charge  coef¬ 
ficients  while  the  self  charge  depends  on  the  absolute  values.  While  the 
vectors  of  the  inverted  matrix  composed  of  the  above  equations  are  linearly 
independent  the  individual  coefficients  are  now  statistically  dependent. 

This  dependence  will  affect  the  error  analysis  seen  in  section  VIII. 


VII.  DETERMINATION  OF  ENHANCEMENT  COEFFICIENTS 


The  determination  of  the  enhancement  coefficients  can  logically  be 
divided  into  separate  determinations  of  field  and  self  charge  coefficients. 
Furthermore  each  set  contains  one  arbitrarily  chosen  absolute  coefficient  (we 
choose  Py  and  Pq)  to  which  the  others  can  be  expressed  as  ratios. 


A.  Self  charge  coefficients 


Once  a  sensitive  enough  range  is  attained,  then  the  ratio  of  the  self 
charge  coefficients  can  be  determined  through  changes  in  the  self  charge 
brought  about  by  varying  the  engine  characteristics  (e.g.  -  varying  rpm  or 
power)  or  by  traversing  small  non-electrif led  cloud  tufts  which  charge  the 
aircraft  by  waterdrop  impaction.  For  the  CV-580  the  charging  by  water  impac¬ 
tion  was  opposite  in  sign  compared  to  the  engine  charging.  An  example  of 
charge  variation  is  given  in  figure  9  where  the  engine  power  was  varied 
sequentially  from  cruise  power  to  low  power  (600  HP)  to  high  power  (2500  HP). 
Apparently  cruise  power  and  high  power  produce  approximately  the  same  char¬ 
ging  rate  on  the  CV-580.  This  may  be  caused  by  the  corona  wires  limiting  the 
maximum  aircraft  potential.  The  mixing  techniques  as  discussed  in  section 
VII. C  can  be  used  as  an  alternative  method  of  determining  these  ratios. 

To  determine  the  absolute  self  charge  coefficient  Pq,  comparisons  must  be 
made  between  in  flight  and  ground  measurements.  The  effects  of  image  charges 
as  the  aircraft  approaches  the  ground  cause  artificial  enhancement  of  the  Pq 
coefficient.  Ground  measurements  of  the  fields  observed  upon  the  application 
of  a  known  DC  voltage  to  the  aircraft,  while  the  wheels  were  resting  on  bases 
of  1/A  inch  (6.35  mm)  plywaod  and  1/16  inch  (1.59  mm)  electrically  insulating 
polyethylene  mats,  were  performed  at  Patuxent  Naval  Air  Station.  The  ratio 
of  the  measured  field  to  the  known  DC  voltage  relates  aircraft  potential  to 
the  field  in  the  presence  of  the  image  aircraft.  The  effect  of  the  image  can 
be  evaluated  from  examination  of  data  recordings  taken  during  the  landing 
process.  Immediately  prior  to  touchdown,  the  image  is  identical  to  that  seen 
on  the  insulating  mats.  Since  the  self  charge  was  dominant  during  the 
landing  approaches,  field  variations  can  be  ignored.  Thus  Pv,  in  terms  of 
the  equivalent  voltage  produced  by  Q,  can  be  defined  for  the  PORT  (p)  mill 
as : 


Pv  =  (Epa/Va)  =  (Epa/Epl)*(Epl/Vl)*(Vl/Va)  (9) 


where  Epa,  Epl,  Va,  and  VI  are  the  fields  at  the  mill  and  potentials  far 
above  the  ground  (where  image  charges  are  negligible)  at  altitude  (a)  and 
just  before  landing  (1).  Assuming  that  the  self  charge  Q  is  conserved  during 


a  landing  approach  (fig.  10),  then  not  only  is  (Epa/Epl)  a  valid  measurement, 
but  we  have  the  following  additional  result: 

Q  =  Ca  Va  =  Cl  VI  (10) 

where  Ca  and  Cl  are  the  capacitances  in  air  and  at  landing.  The  Patuxent 
ground  measurements  on  the  insulated  mats  yield  (Epg/Vg)  which  is  assumed  to 
be  identical  to  (Epl/Vl),  where  Epg  is  the  field  at  the  mill  produced  by  the 
DC  voltage  Vg.  Thus  Pv  becomes: 

Pv  =  (Epa/Epl) *( Epg/Vg )*(Ca/Cl)  (11) 

Note  that  Pv  instead  of  Pq  has  been  found  (Pq  Q  and  Pv  V  are  equivalent). 

The  capacitance  ratio  Ca/Cl  can  be  determined  by  modeling  of  the  air¬ 
craft.  An  expression  for  the  capacitance  of  a  conducting  sphere  above  an 
infinite  conducting  plane  has  been  derived  [Morse  +  Feshbach,  1953] .  It  is 
obvious  that  the  sphere  underestimates  the  true  capacitance  of  an  aircraft  at 
small  heights.  Another  simple  model  which  should  overestimate  the  capaci¬ 
tance  of  the  aircraft  at  small  heights  is  a  circular  disc  parallel  to  a 
ground  plane  [Andersen,  1983].  At  small  heights,  it  is  possible  to  use  the 
parallel  plate  formula  for  capacitance  with  little  error;  while  the  capaci¬ 
tance  of  the  disk  at  infinite  height  is  available  in  analytic  form  [Terman, 
1943].  The  capacitance  variation  with  height  for  the  sphere  and  disk  models 
is  plotted  in  fig.  11  normalized  to  the  values  at  infinite  height.  Both 
sphere  and  disk  were  assigned  a  radius  of  6.76  meters  which  gives  the  same 
horizontal  cross  sectional  area  as  the  CV-580  aircraft.  The  aircraft  height 
above  the  ground  at  the  lowest  point  on  the  bottom  of  the  fuselage  of  1.54 
meters  is  indicated  in  figure  11  by  the  arrow  and  nearly  coincides  with  the 
intersection  point  of  the  two  models.  The  capacitance  ratio  at  this  inter¬ 
section  point  is  0.54. 

As  noted  earlier,  a  value  for  Pv  has  been  obtained.  Pq  is  available  if 
an  absolute  aircraft  capacitance  is  known.  Pv,  which  relates  self  charge 
field  to  voltage  on  the  aircraft,  is  certainly  a  useful  quantity;  but  Pq 
could  be  estimated  from  the  752  pF  capacitance  at  infinity  of  a  sphere  of 


CAPACITANCE 


Normalized  capacitance  values  for  conducting  sphere  and  circular 
disk  as  functions  of  height.  The  height  (1.54  m)  of  the  lowest 
point  on  the  underside  of  the  CV-580  fuselage  is  indicated. 


6.76  meter  radius  (this  estimate  is  probably  high  -  the  disk  model  gives  a 
value  of  479  pF) .  The  measured  values  for  each  factor  in  equation  11  are: 


which  yields 


Epg/Vg  =  2.044  (V/m)/V 
Epa/Epl  =  1.273 

Ca/Cl  =  0.54 

Pv  =  1.41  (V/m) /V 


Corresponding  values  for  Pq  (=  Pv/Ca)  would  be  1875  (V/m)//iC  (C=752)  or  2944 
(V/m)//<C  (C=479)  where  a  realistic  value  might  be  the  mean  value  of  2410 
(V/m)//iC.  The  overall  error  of  Pv  is  determined  from  the  error  of  each 
factor  above.  Epg/Vg  is  known  accurately  (<  3Z)  while  the  error  on  the  capa¬ 
citance  is  difficult  to  ascertain.  The  capacitance  ratio  Ca/Cl  is  fairly 
well  known.  This  is  why  we  calculate  Pv  instead  of  Pq.  The  biggest  error 
term  will  probably  be  the  measurement  of  Epa/Epl. 

Note  that  a  distinction  should  be  made  between  Cl  (just  before  touchdown) 
and  Cg  (sitting  on  the  ground).  Cg  was  also  measured  over  a  concrete  ramp  at 
Patuxent  with  a  value  of  3010  pF.  Given  the  capacitance  ratio  of  Ca/Cl  = 

0.54  and  the  estimates  from  the  models  for  Ca  then  it  is  found  that  Cl  has 
the  limits  of  887  to  1393  pF  which  is  a  factor  of  two  to  three  different  from 
Cg.  The  difference  is  probably  due  to  the  distortion  of  the  tires  when  sit¬ 
ting  on  the  ground  which  lowers  the  aircraft  and  provides  a  larger  footprint 
of  the  wheels.  Both  effects  tend  to  increase  the  measured  capacitance. 

B.  Absolute  field  coefficient  Py 

With  the  instrumentation  on  this  aircraft  and  because  of  the  .igh  self 
charge  it  must  be  pointed  out  that  the  fair  weather  field  can  only  be  meas¬ 
ured  by  difference  techniques  (either  the  field  must  change  quickly  or  air¬ 
craft  maneuvers  must  be  performed);  therefore,  comparisons  with  a  ground  site 
were  not  practical.  This  leaves  two  options  of  absolute  calibration:  1) 
direct  comparison  with  an  aircraft  calibrated  to  measure  the  vertical  field 
or  2)  the  profile  technique  where  a  known  value  is  assumed  for  the 


ionospheric  potential.  The  profiling  technique  is  defined  by  measuring  the 
field  at  regular  altitudes  from  near  ground  height  to  the  ceiling  limits  of 
the  aircraft.  These  measurements  can  then  be  approximately  integrated  to 
obtain  the  ionospheric  potential.  A  direct  comparison  with  another  aircraft 
requires  close  formation  flying.  It  can  eliminate  the  need  for  quickly  cov¬ 
ering  multiple  altitudes  but  has  the  distinct  disadvantages  of  relying  on 
another  aircraft  calibration  and  on  the  coordination  of  aircraft  positions 
and  communication  between  two  research  groups.  The  profile  technique  must  be 
done  quickly  and  depends  critically  on  the  weather  conditions. 

During  direct  comparisons ,  Airborne  Research  Associates’  Baron  aircraft, 
instrumented  with  four  radioactive  potential  probes  (two  on  each  wing  sepa¬ 
rated  vertically  by  1.45  m) ,  performed  joint  maneuvers  with  the  CV-580.  The 
Baron's  calibration  is  based  on  low  level  passes  over  an  instrumented  ground 
site.  Its  accuracy  is  probably  between  5  and  10  percent.  Since  this  is 
lower  than  our  calculated  Py  uncertainty  it  cannot  be  a  major  factor  in  our 
overall  uncertainty.  Also,  spatial  errors  in  the  fields,  due  to  the  air¬ 
crafts  occupying  different  airspaces,  are  probably  much  larger.  The 
Baron’s  contribution  was  essentially  to  measure  the  vertical  fair  weather 
field  while  the  CV-580  did  rolls. 

Given  a  vertical  fair  weather  field  E,  during  a  roll  maneuver  0  (positive 
for  right  banks  -  see  fig.  4)  an  Ey  field  appears  from  zero  ass 

Ey  =  E  sin(0) 

therefore  Py  can  be  calculated  from  left  and  right  banks  as 

Y(fli)  -  Y (02) 

Py  - -  (12) 

E  (sin(#i)  -  sin(#2)) 

Normally  6 j  and  6 £  are  nominally  either  *  30  or  *  45  degrees  respectively. 
Since  the  angle  of  roll  varies  somewhat  during  each  maneuver,  it  is  probably 
better  to  divide  each  point  within  the  roll  by  its  appropriate  angle,  aver¬ 
age,  and  then  difference  with  the  values  in  the  other  angle.  E  is  obtained 
easily  from  the  other  comparison  aircraft. 


For  the  profile  technique,  an  uncalibrated  measurement  of  the  vertical 
field  must  be  integrated  with  respect  to  height  h  to  the  ionsphere  yielding  a 
value  for  the  ionospheric  potential.  Py  can  be  calculated  from: 

/  [Y(51)  -  Y(02)]  dh 

Py  =  -  (13) 

V  (sin(0!)  -  sin(02)) 

where  the  ionospheric  potential  V  has  been  previously  measured.  Since  the 
aircraft  has  limits  on  both  the  lowest  and  highest  altitudes  of  flight,  a 
typical  profile  (fig.  12)  must  be  extrapolated  to  both  the  ground  and  the 
ionosphere.  The  ground  field  was  assumed  to  be  the  same  as  the  lowest  meas¬ 
urement  (200  ft.)  while  the  ionosphere  value  is  corrected  by  the  slope  above 
the  inversion  level  (exponential  tail).  The  apparent  vertical  potential 
plotted  on  top  in  figure  12  is  the  quantity  Py  V.  The  true  ionospheric 
potential  was  calculated  from  a  compilation  of  yearly  ionospheric  potential 
(fig.  13)  by  Markson  [1985].  Eliminating  the  three  outlier  points  around 
year  1960,  the  mean  is  2A6  KV  *  7X.  The  value  of  Py  V  must  be  corrected  for 
the  well  known  diurnal  variations  as  expressed  by  the  Carnegie  curve  (fig. 

14)  [Ault  and  Mauchly,  1926]  [Whipple,  1929].  Four  profiles  were  accom¬ 
plished  during  the  summers  of  1983  and  198A  over  the  ocean  off  the  coasts  of 
Cape  Canaveral,  Florida  and  Atlantic  City,  N.  J. 

Systematic  differences  might  exist  between  the  direct  comparison  and  the 
profile  technique.  The  low  altitude  calibration  of  the  Baron  aircraft  could 
have  possible  problems  of  image  charge  effects.  Furthermore,  altitude  and 
airflow  variations  might  cause  distortion  of  the  ion  plumes  from  the  radio¬ 
active  potential  probes  thus  affecting  the  enhancement  coefficient.  The 
profile  technique  has  uncertainties  in  the  assumed  value  of  the  ionospheric 
potential,  problems  with  accepting  diurnal  variations  on  a  daily  basis,  and 
errors  in  the  extrapolations  to  both  the  ground  and  the  ionosphere.  Problems 
common  to  both  techniques  are  errors  in  reading  the  data  recordings,  uncer¬ 
tainties  in  the  altitude  and  attitude  of  the  aircraft,  nonstationar ity  of  the 
fair  weather  field,  and  spatial  inhomogeneities.  In  figure  15,  where  the 
profile  and  direct  comparison  are  compared,  a  histogram  of  Py  shows  no  signi¬ 
ficant  differences.  Py  was  calculated  from  this  data  set  as  26.5  *  1A.2%. 
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DISTRIBUTION  OF  Py  COEFFICIENT 


Histogram  of  absolute  Py  coefficient.  The  profile  and  direct 
comparison  methods  are  indicated  (see  text). 


C.  Field  ratio  coefficients 


This  section  deals  with  calculating  the  other  field  coefficients  from  the 
knowledge  of  Py.  It  should  be  pointed  out  that  for  determining  the  ratio  of, 
for  example,  Py/Fz,  the  absolute  fair  weather  field  E  need  not  be  known.  In 
this  section,  however,  it  is  assumed  that  Py  is  already  known,  thus,  E  can  be 
determined  from  roll  but  not  pitch  maneuvers.  In  calculating  the  ratios, 
some  error  causes  will  be  eliminated.  There  will  be  no  errors  with  respect 
to  diurnal  corrections,  value  of  the  ionospheric  potential,  and  extrapola¬ 
tions.  Remaining  errors  are  due  to  non-stationarity  and  inhomogeneity  of  the 
fair  weather  field,  uncertainties  in  attitude,  changes  in  altitude  during 
pitch  maneuvers,  and  statistical  fluctuations  in  the  data  (mostly  due  to  self 
charge  variations). 

Due  to  the  symmetry  in  the  wingtips,  the  measurement  of  the  derived  Ey 
field  is  uncontaminated  by  self  charge  variations.  For  the  belly  mills  and 
the  derived  quantity  Q' ,  techniques  must  be  devised  to  eliminate  the  more 
severe  contamination  from  self  charge  variations.  All  the  field  coefficients 
except  Py  are  affected.  A  high  correlation  exists  between  Q' ,  mostly  a 
measure  of  self  charge,  and  both  of  the  belly  mills.  The  mixing  of  the  self 
charge  as  measured  by  Q'  into  the  belly  mills  takes  advantage  of  this  corre¬ 
lation  by  largely  eliminating  self  charge  from  the  belly  mills.  By  choosing 
mixing  parameters  a  and  such  that  known  self  charge  variations  vanish,  the 
belly  mills  are  redefined  as: 


Ef '  5  Ef  -  a  Q' 

-  (Fx 

-  a  Px)Ex  +  (Fz 

-  a  Pz)Ez 

(14) 

Et'  5  Et  -  p  Q' 

*  (Tx 

-  fi  Px)Ex  +  (Tz 

-  P  Pz)Ez 

(15) 

where  the  contaminations  due  to  Px,  and  Pz  are  left  to  be  estimated  (it  will 
be  shown  later  that  Px  and  Pz  are  close  to  zero).  The  values  for  a  and  /) 
were  0.149  and  0.258  respectively.  As  mentioned  earlier,  these  values  can  be 
used  to  calculate  the  ratio  of  the  self  charge  coefficients. 

Each  of  the  belly  mill  enhancement  coefficients  can  be  determined  from  a 
specific  maneuver.  Fz  and  Tz  are  determined  from  roll  maneuvers  while  Fx  and 
Tx  are  determined  from  pitch  maneuvers.  Since  pitch  maneuvers  cannot  measure 


E,  *  set  of  roll-pitch  maneuvers  oust  be  performed  sequentially.  The  normal 
procedures  were  s  period  of  straight  and  level  flight  followed  by  consecutive 
45  degree  left  and  right  banks,  another  period  of  straight  and  level,  a 
series  of  consecutive  climb-dives  (*  15  degree),  a  third  period  of  straight 
and  level,  consecutive  45  degree  left  and  right  banks,  and  a  final  period  of 
straight  and  level.  This  allows  E  to  be  calculated  during  the  rolls  for  use 
in  the  climb-dive  maneuvers.  With  the  above  expressions  for  the  belly  mills 
uncontaminated  by  self  charge,  the  four  coefficients  Fx,  Fz,  Tx,  and  Tz  can 
be  determined  as: 

Ef'(fi)  -  Ef'(O) 

Fz  -  a  Pz  -  -  (16) 

E  (cos 0  -  1) 

Et'(0)  -  Et v (0) 

Tz  -  P  Pz  -  -  (17) 

E  (cos 8  -  1) 

Ef'(^)  -  Ef*(^2) 

Fx  -  a  Px  -  -  (18) 

E  ( sin^i  -  sin^) 

Et ' (^i )  -  Et'(*z) 

Tx  -  p  Px  -  -  (19) 

E  (sin^i  -  sin^2) 

where  f  is  the  pitch  angle.  Note  that  *  15  degree  pitch  angles  are  better 
than  30  degree  dives  since  they  can  be  performed  in  sequences,  less  altitude 
variation  occurs,  the  pilots  are  more  willing  to  perform  them,  and  the  rela¬ 
tive  change  is  greater  (2  sinl5  >  sin30). 

The  mixing  techniques  described  above  are  illustrated  in  fig.  16-19  for 
the  Fx,  Tx,  Fz,  and  Tz  coefficients.  The  improvement  in  the  data  qua  lity  is 
clearly  evident  when  the  top  (Q‘  not  removed)  is  compared  to  the  bottom  plot 
(Q'  removed).  The  data  traces,  at  least  for  Fx,  Tx,  and  Fz  (fig.  16-18)  do 
exactly  what  is  expected.  Roll  maneuvers  for  the  belly  mills  will  produce  a 
reduction  in  the  voltage  seen  by  the  mill  in  both  a  left  and  right  bank  (due 
to  the  evenness  of  the  cosine  function).  Pitch  maneuvers,  on  the  other  hand, 
will  produce  a  sinusoidal  behavior  depending  on  whether  a  climb  or  a  dive  is 
being  performed.  The  changes  seen  in  the  Fx  and  Tx  traces  should  be  of  oppo¬ 
site  polarity  while  the  changes  in  the  Fz  and  Tz  traces  should  be  of  same 


RELATIVE  SYSTEM  OUTPUT  (V)  RELATIVE  SYSTEM  OUTPUT  (V) 

FORWARD  -  WITH  Q 'MIXING  FORWARD  -  NO  Q 'MIXING 


CALIBRATION  FLIGHT  (08/12/85) 


TIME  (sec) 


-5r  2Q1 


TIML  (sec) 


Fig.  16  - 


Fx  coefficient  -  mixing  of  Q'  into  FORWARD  mill  during  climb-dives 
to  minimize  Q  variations  (top  -  no  mixing,  bottom  -  with  mixing). 
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7  -  Tx  coefficient  -  mixing  of  Q'  into  TAIL  mill  during  climb-dives  to 
minimize  Q  variations  (top  -  no  mixing,  bottom  -  with  mixing). 
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Fig.  18  -  Fz  coefficient  -  mixing  of  Q'  into  FORWARD  mill  during  rolls  to 
minimize  Q  variations  (top  -  no  mixing,  bottom  -  with  mixing). 
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Fig.  19  -  Tz  coefficient  -  mixing  of  Q'  into  TAIL  mill  during  rolls  to 

minimize  Q  variations  (top  -  no  mixing,  bottom  -  with  mixing). 
Because  of  the  poorer  data  quality,  the  region  used  for  the 
analysis  is  indicated  (see  text). 
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As  can  be  seen  from  fig.  19,  the  measurement  for  the  Tz  coefficient  is 
not  readily  apparent.  Apparently,  for  this  case,  either  the  self  charge 
removal  was  inadequate  or  the  signal  is  inherently  noisier  due  perhaps  to 
space  charge  plumes  emitted  from  the  engines  drifting  past  the  TAIL  meter. 
Since  good  quality  traces  over  the  time  of  the  roll  maneuver  were  not  easily 
identifiable,  the  data  analyzed  was  the  small  region  (as  indicated)  between 
the  left  and  right  banks  as  the  roll  angle  passes  through  straight  and  level. 
Even  with  this  approach,  many  traces  were  still  unusable,  therefore,  two 
criteria  were  adopted.  Since  Q  variations  are  evidently  very  important,  in 
the  first  criterion,  cases  were  chosen  such  that  no  severe  changes  in  Q' 
occurred  during  the  short  period  of  time  of  the  indicated  measurement  region. 
The  second  criterion,  requiring  equivalent  amplitudes  on  both  sides  of  the 
straight  and  level  after  a  return  to  the  45  degree  angle,  helps  eliminate 
linear  trends. 

Other  than  the  suggestion  of  the  exhaust  plumes,  there  is  no  apparent 
reason  as  to  why  the  measurement  of  Tz  should  be  bad.  There  is  a  factor  of 
approximately  2  reduction  in  the  change  in  signal  seen  in  a  pitch  versus  a 
roll  maneuver  (2  sinl5  vs.  cos45  -  1);  however,  the  same  roll  maneuver  gives 
reasonable  results  for  the  FORWARD  mill.  Also,  the  same  TAIL  mill,  during  a 
pitch  maneuver  gives  accurate  results  for  the  Tx  coefficient.  It  is  also 
reasonable  to  assume  that  the  value  obtained  for  Tz  should  be  close  to  the 
value  obtained  for  Fz.  As  is  seen  later,  this  is  indeed  borne  out.  There¬ 
fore,  even  though  the  data  quality  is  suspect,  the  value  obtained  appears  to 
be  about  what  would  be  expected. 

The  values  of  both  the  Px  and  the  Pz  coefficients  must  be  known  to  elimi¬ 
nate  their  contamination  in  the  above  calculations.  In  order  to  measure 
these  coefficients,  we  must  examine  the  derived  quantity  Q'  (eq.  4),  which 
measures  mostly  self  charge,  during  both  roll  and  pitch  maneuvers.  In  prin¬ 
ciple,  Px  and  Pz  can  be  determined  from  pitch  and  roll  maneuvers  in  a  similar 
manner  as  Fx,  Tx,  Fz,  and  Tz  from: 

Q'(^l)  -  Q'(fo> 

Px  =  - -  (20) 

E  (sin^  -  sin^) 

Q'(0)  -  Q'(0) 

E  (cos 9  -  1) 


Pz 


(21) 


Unfortunately  Q  is  strongly  contaminated.  An  independent  measure  of  the 
self  charge  is  required.  Fortunately,  since  during  fair  weather  there  is 
only  a  field  component  (E)  in  the  vertical  direction,  a  combination  of  the 
FORWARD  and  TAIL  mills,  now  with  some  knowledge  of  the  belly  mill  coeffi¬ 
cients,  gives  the  required  independent  measure.  For  instance  during  a  pitch 
maneuver  the  FORWARD  and  TAIL  mill  responses  are: 

Ef  =  (Fx  sin^  +  Fz  cos^)E  +  Fq  Q 
Et  =  (Tx  sin^  +  Tz  cos^)E  +  Tq  Q 

Solving  these  equations  for  Q  gives: 

(Tx  sin^  +  Tz  cos^)Ef  -  (Fx  sin^  +  Fz  cos^)Et 

q  =  -  (24) 

(Fq  Tx  -  Tq  Fx)sin^  +  (Fq  Tz  -  Tq  Fz)cos^ 

Unfortunately,  since  Tz  and  Fz  are  the  same  sign,  this  equation  for  Q  gives 
inaccurate  results  since  the  denominator  for  small  angles  (FqTz  -  TqFz)  is 
close  to  zero.  The  result  is  the  near  elimination  of  Q  along  with  E  from 
equations  22  and  23,  where  Q  is  what  was  to  be  calculated.  The  above  solu¬ 
tion  essentially  differences  two  mills  that  are  on  the  same  side  of  a  nodal 
curve  (Ez  direction)  where  it  is  difficult  to  separate  the  effect  of  the  self 
charge  and  the  field  contribution.  Better  mill  placement  would  be  similar  to 
the  wingtip  mills,  which  are  symmetrically  on  opposite  sides  of  the  Ey  nodal 
curve.  Since  the  belly  mills  are  on  opposite  sides  of  the  Ex  nodal  curve, 
then  fields  in  only  the  x  direction  can  effectively  be  eliminated  (e.g.  - 
terms  involving  FxEt  and  TxEf).  In  this  case,  the  two  field  coefficients  Fx 
and  Tx  are  now  opposite  in  sign.  The  result  is: 

Tx  Ef  -  Fx  Et  +  (Fx  Tz  -  Tx  Fz)cos^  Ez 

Q  =  -  {25) 

Tx  Fq  -  Fx  Tq 

The  last  term  is  kept  since,  even  though  it  is  small  compared  to  Q,  it  is  of 
the  same  order  of  magnitude  as  the  coefficient  (in  this  case  Px)  being  meas¬ 
ured.  If  the  angles  of  pitch  are  the  same  on  both  sides  (i.e.  -  *  13 
degrees)  then  for  the  pitch  case  the  last  term  will  cancel.  For  roll  maneu¬ 
vers,  a  belly  mill  should  be  differenced  with  a  mill  on  top  of  the  fuselage. 
Since  no  such  mill  is  available,  a  similar  approach  as  with  the  pitch  maneu¬ 
vers  is  taken.  The  result  is  the  same  as  equation  25  except  that  ^  is 
replaced  by  9.  Since,  according  to  eq.  21,  the  comparison  is  made  between 


(22) 

(23) 


banks  and  level  flight,  the  last  term  cannot  be  dropped.  Typical  examples  of 
the  elimination  of  the  self  charge  from  Q'  are  shown  in  figures  20  and  21. 
Each  figure  shows  both  Q'  (top)  and  Q'-PqQ  (bottom),  where  Q  is  found  from 
equation  25.  In  the  data  analysis  further  corrections  have  been  made  to 
compensate  for  linear  trends.  We  can  now  substitute  Q'-Pq  in  place  of  Q' 
into  eqs.  20  and  21  to  find  Px  and  Pz.  Now  that  values  for  Px  and  Pz  are 
known,  their  contamination  in  the  belly  coefficients  can  be  eliminated. 

Some  further  discussion  is  warranted  on  the  trends  in  Q  that  appear 
during  climb-dive  maneuvers  (fig.  20).  They  are  not  due  to  a  real  measure¬ 
ment  of  Px  since  the  variations  fail  to  align  with  the  pitch  angle  (there 
appears  to  be  some  phase  offset).  Also  they  can  easily  be  identified  as  self 
charge  variations  from  the  independent  belly  mill  measurements.  Apparently 
the  self  charge  on  the  aircraft  is  always  corona  limited.  Rates  of  corona 
discharge  are  sensitive  to  airspeed,  pressure,  humidity,  temperature,  and  the 
absolute  magnitude  and  direction  of  the  external  field.  A  simple  model  of 
the  sensitivity  of  the  corona  discharge  for  an  isolated  point  in  space  under 
quasi-static  conditions,  considering  airspeed  only,  has  been  developed 
[Chalmers,  1967] J 

ic  =  A  V  u  (26) 

where  ic,  A,  V,  and  u  are  the  corona  discharge  current,  proportionality  con¬ 
stant,  aircraft  potential,  and  airspeed,  respectively.  Here  the  assumptions 
have  been  made  that  the  engines  maintain  a  constant  charging  current  and  self 
charge  variations  are  caused  only  by  the  corona  wires.  Determining  the  air¬ 
speed  simply  from  gravitational  acceleration  and  ignoring  drag,  a  solution 
for  the  shape  of  V  can  be  found: 

dc 

V  =  -  (27) 

A  (Wq  +  g  /  sin(^(t)  -  fo)  dt> 

where  Wg  =  90.1  m/s  (175  Kts)  is  the  aircraft  cruising  speed,  g  is  the  gravi¬ 
tational  constant,  and  =  +3.5*  is  the  cruising  pitch  angle.  A  simple 
check  on  the  validity  of  the  results  predicted  by  this  formula  can  be  made  by 
examining  the  airspeed.  The  maximum  (157  m/s  =  305  Kts)  and  minimum  (66.5 
m/s  =*  129  Kts),  as  derived  from  the  data,  should  not  greatly  exceed  the  red- 
line  (250  Kts  at  low  altitude)  and  stall  (90  -  95  Kts)  speeds  for  the  air¬ 
craft.  Since  the  drag  term  has  been  ignored  these  values  are  probably 
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Fig.  20  -  Px  coefficient  -  subtracting  Pq  Q  as  measured  independently  from 
the  FORWARD  and  TAIL  mills  during  pitch  maneuvers  (see  text). 
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reasonable.  By  choosing  ic/A  such  that  the  theory  and  data  approximately 
match,  the  results  In  £igure  22  are  obtained.  Ic/A  for  the  three  cases  Is 
the  same  and  was  chosen  to  give  the  best  overall  fit.  Obviously  there  Is 
more  going  on  than  this  simple  model  can  explain;  however,  it  does  suffice  to 
reproduce  the  overall  shape.  It  should  be  noted  that  this  model  was  not  used 
in  determining  the  enhancement  coefficients. 

Since  all  coefficients  are  now  known,  further  refinements  can  be  made  by 
using  the  derived  (from  solving  the  matrix  -  eq.  8)  Q  to  eliminate  the  self 
charge  variation  in  all  channels  except  Q' .  An  iterative  approach  can  be 
taken.  The  values  from  this  technique  taken  along  with  Px  and  Pz  as  dis¬ 
cussed  above  are  the  final  results  for  the  enhancement  coefficients. 

The  relative  accuracy  of  an  enhancement  coefficient  depends  both  on  the 
field  mill  location  and  on  the  type  of  maneuvers  performed.  By  far  the  worst 
coefficients  are  Px,  Pz,  and  Tz.  Among  these  the  accuracy  of  Tz  probably 
matters  the  most.  This  is  true  because  its  relatively  larger  magnitude 
affects  the  inverted  matrix  more.  Also,  from  the  error  analysis  In  section 
VIII,  it  appears  that  changes  in  the  error  of  this  coefficient  affect  the 
errors  in  the  field  more  than  changes  in  the  error  of  either  Px  or  Pz.  In 
general,  pitch  maneuvers  are  more  difficult  to  perform  than  roll  maneuvers. 
The  pitch  angle  is  harder  to  maintain  (it  is  essential  to  record  it),  varia¬ 
tions  in  altitude  affect  the  measured  external  field;  and,  as  alluded  to 
previously,  changes  in  the  aircraft  speed  can  change  the  self  charge.  Fur¬ 
thermore,  a  pitch  maneuver  cannot  attain  the  larger  angles  of  a  roll  maneuver 
(*  15*  versus  *  45*). 

Several  attempts  were  made  to  obtain  corroborative  information  of  the 
enhancement  coefficient  values.  Stable  vertical  high  field  conditions  would 
give  much  cleaner  data  traces  because  the  relative  changes  seen  during  a 
maneuver  would  be  larger  while  the  self  charge  variations  would  presumable 
stay  the  same.  Near  the  anvil  of  a  thunderstorm,  a  high  field  would  presum¬ 
ably  exist.  There  would,  however,  be  questions  as  to  the  stability.  Never¬ 
theless,  anvils  were  tried;  however,  they  proved  mostly  unusable  because 
conditions  varied  to  much.  Data  that  appeared  reasonable  initially  turned 
out  to  be  something  else;  for  instance,  the  deflection  observed  during  a  roll 
maneuver,  upon  further  scrutiny,  failed  to  align  with  the  start  and  the  end 
of  the  roll.  Presumably  the  aircraft  was  not  far  enough  away  from  the  turret 
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of  the  thunderstorm,  causing  a  non-vertical  field.  What  anvil  results  are 
available  do,  however,  usually  appear  to  give  at  least  the  right  sign  of  a 
particular  coefficient. 


Another  approach  would  perhaps  be  an  examination  of  field  changes  from 
lightning  that  does  not  strike  the  aircraft.  Unfortunately,  because  the 
aircraft  is  in  corona  all  the  time,  the  large  field  changes  needed  for  an 
adequate  measurement  also  cause  changes  in  the  self  charge  during  the  light¬ 
ning  discharge.  The  time  constant  for  self  charge  is  estimated  from  a  take¬ 
off  procedure  (fig.  23)  as  less  than  5  seconds.  This  technique  might  work 
for  more  sensitive  ranges  es  we  approach  a  storm  from  farther  away  since  now 
the  fields  are  less  and  might  not  affect  the  self  charge.  Corroborative 
information  might  be  obtained  on  the  Tz  coefficient.  The  problem  of  the 
aircraft  in  corona  was  not  recognized  until  operations  had  been  completed. 

A  similar  technique  using  fast  field  change  radiation  fields  has  been  sug¬ 
gested  by  Beasley  and  Master  [1982].  It  would  require  the  installation  of  a 
fast  field  antenna  at  the  same  location  as  the  field  mill  system.  The  magni¬ 
tude  of  the  field  changes  would  have  to  be  measured  at  a  ground  site.  This 
approach  might  also  work  for  slow  field  changes.  A  new  design  by  Zaepfel 
[1986],  which  combines  a  field  mill  and  slow  antenna  in  one  unit,  could  prob¬ 
ably  suffice  for  the  slow  field  change  measurement.  Finally,  modelling  of 
the  aircraft  geometry  could  be  examined  as  a  comparison  with  the  experiment¬ 
ally  determined  coefficients. 


D.  Values  of  coefficients 


The  final  average  enhancement  coefficients,  along  with  their  sample  stan¬ 
dard  deviations,  are 


Px-  .503*1.79 

Py-  26.5*3.76 

Pz-  0  *  1.95 

Pv» 

1.41*0.07 

Sx«  .503*1.79 

Sy— 26.5*3.76 

Sz-  0  *  1.95 

Sv= 

1.41*0.07 

Fx«  1.69*. 391 

Fz— 4. 24*.  854 

Fv» 

0.25*0.01 

Tx— 3. 24*.  340 

Tz— 3.62*1.25 

Tv- 

0.32*0.02 

(28) 
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Since  the  TAIL  mill  is  further  aft  of  the  wings,  the  presumed  Ex  nodal  curve 
than  the  FORWARD  mill  is  ahead  of  them,  it  is  reasonable  to  expect  that  Tx 
will  be  of  opposite  sign  and  larger  magnitude  than  Fx.  Similarly,  since  the 
TAIL  mill  is  slightly  closer  to  the  Ez  nodal  curve  than  the  FORWARD  mill,  Tz 
should  be  of  same  sign  and  slightly  less  in  magnitude  than  Fz.  The  wingtip 
coefficients  Px,  Sx,  Pz,  and  Sz  should  be  close  to  zero  because  of  the  proxi 
mity  of  the  wingtip  meters  to  the  wingtip  nodal  points.  Also  indicated  is 
the  stronger  sensitivity  of  the  wingtips  to  the  self  charge  than  the  belly 
mills. 

Since  each  of  the  coefficient  values  quoted  in  eq.  28  is  actually  the 
average  of  a  number  of  independent  samples,  a  more  realistic  uncertainty 
t*>uld  be  s/Jn  ,  where  s  is  the  standard  deviation  of  the  individual  samples 
and  N  is  the  number  of  degrees  of  freedom.  It  is  difficult  to  determine  the 
appropriate  value  of  N,  however,  because  of  the  poor  quality  of  some  of  the 
samples  -  especially  those  for  Px,  Pz,  and  Tz.  Therefore,  the  more  conser¬ 
vative  approach  of  using  the  sample  standard  deviations  has  been  chosen  for 
the  analysis  below.  The  error  estimates  from  this  conservative  approach  are 
probably  more  reasonable  since  any  reduction  in  uncertainties  might  be  com¬ 
pensated  by  the  neglect  of  other  errors  such  as  symmetry  assumptions. 

VIII.  ERROR  ANALYSIS 

Determination  of  the  errors  in  the  field  components  and  the  aircraft 
potential  is  not  trivial  because  of  the  non-linear -matrix  inversion  process. 
In  eq.  2,  it  was  assumed  that  the  matrix  elements  are  statistically  indepen¬ 
dent.  As  can  be  seen  by  the  explicit  dependence  on  the  enhancement  coeffi¬ 
cients  in  the  inverted  matrix  (eqs.  5-8),  the  inverted  coefficients  are  now 
statistically  dependent,  thus,  cross  terms  must  be  considered.  In  order  to 
address  the  problem,  a  Monte  Carlo  method  has  been  used.  First,  it  must  be 
noted  that  due  to  the  symmetry  assumptions  for  the  wingtip  mills,  Py  enters 
only  into  the  uncertainty  associated  with  an  Ey  field.  The  uncertainties 
resulting  from  the  method  will  depend  on  the  assumed  external  field 
conditions . 


Th«  details  of  the  analysis  method  entail  assuming  an  external  field  and 
aircraft  potential.  The  raw  fields  at  the  mill  (Ep,  Es,  Ef,  and  Et)  can  then 
be  calculated  from  the  mean  value  of  the  enhancement  coefficients  given  in 
eq.  28.  One  then  constructs  Gaussian  generators  with  which  the  experiment¬ 
ally  detarmined  enhancement  coefficients  can  be  varied  over  their  one  sigma 
range  (the  standard  deviation  in  eq.  28).  There  will  be  one  such  generator 
for  each  of  the  enhancement  coefficients.  A  perturbed  matrix,  derived  from 
these  generators,  is  inverted  to  obtain  a  set  of  perturbed  inverted  coeffi¬ 
cients.  From  this  perturbed  matrix  and  from  the  previously  calculated  fields 
at  the  mill,  the  perturbed  external  fields  and  aircraft  potential  can  be 
calculated.  The  method  was  repeated  1000  times  until  the  statistics  were 
stable.  The  means  and  standard  deviations  for  the  inverted  coefficients,  the 
field  components ,  and  aircraft  potential  were  then  calculated.  Several  sets 
of  input  field  conditions  were  tested  so  that  regions  of  data  validity  could 
be  estimated. 

During  each  inversion,  a  check  was  made  to  insure  that  singularities  did 
not  occur  (determinant  of  zero).  Furthermore,  it  was  determined  where  pos¬ 
sible  singularities  could  occur.  This  was  done  by  holding  all  input  coeffi¬ 
cients  except  one  constant  and  finding  for  what  value  of  the  free  coefficient 
a  singularity  occurred.  These  values  are  given  in  table  1,  where  columns  2 
and  3  are  Identical  to  the  values  in  equation  28.  The  fourth  column  is  the 
value  for  the  particular  coefficient  which  yields  zero  for  the  determinant  of 
matrices  of  the  form  of  equation  28.  The  last  column  is  the  absolute  value 
of  the  difference  between  the  fourth  and  second  columns  divided  by  the  third 
column.  This  ratio  indicates  how  far  out  on  the  wing  of  the  distribution  the 
singularity  occurs.  Assuming  a  normal  distribution,  the  probability  of  the 
worst  case  (7.00)  occurring  is  2.56  X  10"^. 

The  inverted  coefficients  and  their  standard  deviations  are 


Ex 

Ey 

Ez 

V 


Ep  Es  Ef  Et 

(  .008*. 008)  (  .008*. 008)  (  .181*. 060)  (-.212*. 038) 

(  .019*. 003)  (-.019*. 003) 

(  .024*. 006)  (  .024*. 006)  (-.168*. 044)  (-.080*. 033) 
(  .352*. 043)  (  .352*. 043)  (-.065*. 352)  (  .076*. 307) 


(29) 
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TABLE  1  - 

Possible  singularities 
matrix  (See  text) 

for  inversion 

of  enhancement 

coefficient 

ENHANCEMENT 

COEFFICIENT 

COEFFICIENT 

VALUE  FOR 

NUMBER  OF 

COEFFICIENT 

VALUE 

ERROR  (la  ) 

DET  =  0 

a  AWAY 

Px 

0.503 

1.79 

-62.0 

34.9 

Py 

26.5 

3.76 

0 

7.05 

Pz 

0 

1.95 

-20.9 

10.7 

Pv 

1.41 

0.07 

-.011 

20.3 

Fx 

1.69 

.391 

-3.84 

14.1 

Fz 

-4.24 

.854 

1.73 

7.00 

Fv 

0.25 

0.01 

15.8 

1550 

Tx 

-3.24 

0.34 

1.48 

13.9 

Tz 

-3.62 

1.25 

8.88 

10.0 

Tv 

0.32 

0.02 

-12.9 

662 

Note  that  the  errors  of  interest  are  in  the  fields  and  self  charge  which 
depend  not  only  on  the  errors  in  the  matrix  elements  but  also  on  the  magni¬ 
tude  and  direction  of  the  field  and  on  the  self  charge.  The  direction  angles 
are  the  standard  azimuth-elevation  coordinates  as  indicated  in  figure  24. 
(Note  that  these  angles  are  not  the  same  as  in  figure  4).  The  errors  are 
listed  in  table  2  where  some  input  conditions  have  .been  chosen  and  the  rela¬ 
tive  errors  calculated.  The  values  in  parenthesis  are  in  percent  for  the 
errors  of  concern  for  a  particular  case.  For  instance,  if  Ex  is  dominant, 
then  the  relative  errors  in  Ey  and  Ez  do  not  matter.  An  estimate  of  the 
relative  magnitude  of  the  self  charge  contribution  at  each  mill  with  respect 
to  the  field  contributions  (Vmill/Emill)  has  been  determined  as  |V|  times  the 
average  of  the  absolute  values  of  the  four  self  charge  coefficients  divided 
by  |E|  times  the  average  of  the  absolute  value  of  the  twelve  field  coeffi¬ 
cients.  The  result  is 


Vmill 


|V| 

0.152  - 

|E| 


Emill 


(30) 


Azimuth-elevation  coordinate  system  for  field  magnitude  and 
directions. 


rrors  in  the  external  fields  and  aircraft  self  charge  for  some  typical  cases  (both  compone 
irection  are  given) . 
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This  term  does  not  indicate  changes  in  errors  due  to  different  field  direc¬ 
tions.  A  value  of  one  for  this  parameter  should  approximately  indicate  equal 
contributions  of  the  aircraft  potential  and  external  fields  to  the  meters. 
This  case  is  the  fifth  row  where  the  self  charge  is  11. A  times  any  of  the 
field  components.  Our  range  of  operation  for  (Vmill/Emill)  during  non-strike 
conditions  is  approximately  0  to  2.  Not  surprisingly,  it  varies  dramatically 
during  times  of  lightning  attachment. 

Errors  typically  are  about  20%  for  the  magnitude  and  direction  of  the 
external  field  and  for  the  self  charge.  On  first  examination  it  would  appear 
that  the  Ex  field  has  the  highest  error  (row  1);  however,  in  comparing 
approximately  equal  field  contributions  in  all  directions  versus  a  field  in 
mostly  one  direction  (rows  2  and  4),  we  see  that  the  Ex  field  error  lowers 
drastically  (37 %  to  13%)  while  the  Ez  field  error  only  lowers  from  27%  to 
23%.  This  may  be  telling  us  that  the  worst  error  is  really  in  the  Ez  direc¬ 
tion.  This  agrees  with  the  fact  that  both  belly  mills  tend  to  measure  the 
same  sign  of  the  Ez  field  but  opposite  signs  of  the  Ex  field  (opposite  sides 
of  a  nodal  curve  -  thus  difference  techniques  can  be  employed).  It  is  diffi¬ 
cult  to  tell  at  this  point  where  a  limit  should  be  made  in  errors  so  that  one 
component  cannot  be  considered  useful.  This  limit  will  depend  to  some  extent 
on  what  is  required  from  the  data  set  -  in  our  case  we  do  not  yet  know  how 
accurate  we  need  to  be  to  evaluate  the  triggering  hypothesis.  An  attempt  was 
nevertheless  made  to  get  an  indication  of  this  limit. 

In  figure  25,  we  have  plotted  the  error  budgets,  in  terms  of  percent,  for 
Ex,  Ez,  6,  |E|,  and  V,  against  the  (Vmill/Emill)  parameter  for  field  com¬ 

ponents  that  are  identical  in  all  directions.  The  curves  are  plotted  for 
constant  Ex=Ey=Ez=l  while  varying  V.  For  fields  in  other  directions,  the 
error  curves  are  of  approximately  the  same  shape  but  will  be  offset  vertic¬ 
ally  (usually  to  a  lower  error).  The  error  in  Ey  has  not  been  included  since 
it  has  a  constant  value  of  16%.  It  is  seen  that  the  error  on  the  self  charge 
V  is  less  than  30%  for  Vmill/Emill  >  0.5.  The  field  error  terms  appear  to 
start  to  increase  more  rapidly  at  about  Vmill/Emill  =  2.0;  therefore,  all 
components  will  be  considered  valid  if  0.5  £  Vmill/Emill  £  2.0.  This  does 
not  mean  that  the  field  components  are  not  valid  below  0.5  or  the  self  charge 
above  2.0.  Indeed,  this  is  where  these  are  most  valid.  Of  particular  impor¬ 
tance  are  the  magnitude  and  angle  errors  which  are  less  than  25%  in  the  above 
Vmill/Emill  range. 


ERROR  PLOTS  FOR  DERIVED  FIELDS  RND  SELF  CHRRGE 
FOR  R  TYPICRL  EXTERNRL  FIELD 


IX.  IMPROVEMENTS 


The  analysis  of  this  data  has  suggested  some  improvements  for  future 
reference.  It  appears  valuable  to  have  some  redundancy  in  the  number  of 
field  mills  to  add  accuracy  and  reliability  to  the  results.  A  nose  location, 
although  ideal  for  the  Ex  field  component  may  have  severe  problems  with 
impaction  and  the  need  for  a  more  powerful  motor  to  overcome  the  increased 
drag  on  the  rotor  caused  by  the  aircraft  movement.  A  better  location  for  an 
extra  field  meter  might  be  on  the  top  of  the  fuselage  in  order  to  provide  a 
meter  on  the  opposite  side  of  the  Ez  nodal  curve  so  that  the  belly  mills  can 
be  differenced  against  it.  Some  provision  would  have  to  be  made  to  prevent 
this  field  meter  from  filling  with  water.  For  the  CV-580,  since  all  the 
nodal  points  are  not  accessible,  an  optimum  installation  might  include  the 
existing  locations  plus  two  extra  field  meters  on  the  top  of  the  fuselage 
symmetrically  above  the  belly  meters.  This  would  provide  dual  measurements 
of  the  least  accurately  known  fields  Ex  and  Ez.  An  intriguing  possibility  is 
to  combine  a  nose  mounted  cylindrical  mill  with  shutter  mills.  Finally,  an 
active  discharger  (see  Pelton  et.  al.  [1953]),  if  not  too  noisy,  (or  another 
aircraft  with  a  lower  self  charge)  might  aid  in  reducing  the  deleterious  self 
charge  effects  during  calibrations. 


X.  CONCLUSIONS 

The  aim  of  this  work  has  been  the  full  vector  calibration  of  an  aircraft 
instrumented  with  shutter  field  mills.  The  details  involve  experimentally 
determining  the  enhancement  coefficients,  which  when  applied  to  the  fields  as 
seen  at  the  field  meter  locations,  yield  the  true  external  field  and  aircraft 
self  charge.  The  two  most  important  aspects  of  a  successful  calibration  are 
the  initial  field  meter  design  and  the  choice  of  their  locations  on  the  air¬ 
craft.  This  is  probably  not  too  surprising,  but  cannot  be  overemphasized. 

Our  choice  of  calibration  procedure  was  an  experimental  calibration  based  on 
aircraft  maneuvers  in  the  vertical  fair  weather  field;  thus  our  sensitivity 
had  to  be  on  the  order  of  tens  of  volts  per  meter.  With  a  distance  between 
the  rotor  and  stator  of  the  order  of  mm,  contact  potentials  can  be  a  problem. 
A  contact  potential  of  0.1  volts  gives  offsets  in  the  field  of  100  V/m  for  a 
spacing  of  1  mm.  A  previous  study  (Willett  +  Bailey,  1983]  indicates  that 
only  stainless  steel  should  be  used  in  the  sensing  area  in  order  to  minimize 


contact  potentials.  Field  meter  locations  should  be  at  ’nodal  points*  or  at 
the  least  should  be  on  opposite  sides  o£  nodal  curves  (see  section  IV). 
Un£ortunately,  optimum  locations  are  not  always  accessible.  The  third  major 
problem  to  overcome  is  the  deleterious  effects  of  self  charge.  This  problem 
entails  not  only  the  reduction  of  the  absolute  level  but  also  of  the  relative 
fluctuations  in  the  self  charge.  The  relatively  high  magnitude  of  the  self 
charge  seen  aboard  the  CV-580  was  reduced  by  attaching,  in  addition  to  the 
standard  discharge  wicks,  smaller  diameter  (0.003  in)  stainless  steel  wires 
to  the  aircraft.  The  remaining  fluctuations  were  minimized  by  mixing  inde¬ 
pendent  measures  of  the  self  charge  appropriately  into  the  contaminated  field 
meter  outputs. 

Calibration  of  the  CV-580  aircraft  has  been  obtained.  Determination  of 
the  enhancement  coefficients  was  not  trivial  since  the  determination  required 
stability  (or  at  least  a  measurement  of  the  instabilities)  of  both  the  exter¬ 
nal  field  and  the  self  charge  during  a  series  of  maneuvers.  Because  the  self 
charge  fluctuations  mask  the  fair  weather  signal,  maneuvers  are  required  in 
order  to  provide  measurable  differences  in  the  raw  outputs  of  the  mills.  The 
relative  ratio  of  the  enhancement  coefficients  was  determined  without  know¬ 
ledge  of  the  absolute  value.  The  absolute  value  was  then  obtained  by  two 
independent  methods  -  direct  comparison  with  another  aircraft  calibrated  with 
respect  to  the  vertical  field  and  the  profile  technique,  which  involves  inte¬ 
grating  the  field  to  the  ionosphere.  The  absolute  self  charge  coefficients 
were  determined  from  ground  measurements  on  insulating  mats  in  combination 
with  estimates  of  image  effects  from  landing  approaches. 

After  the  enhancement  coefficients  have  been  found,  inversion  of  the 
matrix  will  yield  the  vector  field  components.  These  inverted  matrix  coeffi¬ 
cients  are  given  in  equation  29.  It  should  be  noted  that  the  absolute  value 
of  the  field  components  can  be  determined  with  a  knowledge  of  only  the  ratios 
of  the  self  charge  coefficients.  The  external  field  and  self  charge  errors 
are  estimated  from  a  Monte  Carlo  technique.  The  error  on  the  magnitude  of 
the  external  field  is  less  than  25X  for  an  aircraft  potential  that  is  not 
dominant  (|V|/|E|  <  13.1)  and  the  errors  in  the  direction  angles  are  normally 
less  than  10*.  The  error  in  the  self  charge  is  normally  less  than  20%.  It 
is  suggested  that  if  you  want  to  measure  both  the  field  and  the  self  charge 
at  the  same  instant  with  reasonable  success  then  the  condition  3.28  £  ( V | / ( E | 


£  13.1  must  be  satisfied.  In  a  companion  report  (Anderson  +  Bailey  [1987],  V 
and  | E |  can  be  found  in  figures  66-95.  Thus  |V|/|E|  can  easily  be 
determined. 
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